Mechanical modulation of single-electron tunneling through molecular-assembled 

metallic nanoparticles 
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We present a microscopic study of single-electron tunneling in nanomechanical double-barrier 
tunneling junctions formed using a vibrating scanning nanoprobe and a metallic nanoparticle con- 
nected to a metallic substrate through a molecular bridge. We analyze the motion of single elec- 
trons on and off the nanoparticle through the tunneling current, the displacement current and the 
charging-induced electrostatic force on the vibrating nanoprobe. We demonstrate the mechanical 
single-electron turnstile effect by applying the theory to a gold nanoparticle connected to the gold 
substrate through alkane dithiol molecular bridge and probed by a vibrating platinum tip. 

PACS numbers: 85.35.Gv,85.85.+j,85.65.+h,73.63.-b 



Introduction. — Nanoelectromechanical devices that 
combine mechanics with electronics are of great inter- 
ests for applications in electronics, precision measure- 
ment and sensors*- 1 Among the experimental implementa- 
tions, mechanical single-electron devices that explore the 
interplay between the macroscopic motion of a nanome- 
chanical element and the quantized single-electron trans- 
fer have attracted much attention for revealing new mech- 
anisms of electron transport and schemes of mechanical 
detection at the quantum-limit^ In particular, the shut- 
tle mechanism of quantized charge transfer has been pro- 
posed theoretically 3 ^ and studied experimentally^ utiliz- 
ing a variety of nanomechanical elements. 

The model system considered in the original proposal 3 
for a mechanical single-electron shuttle consists of a small 
metal cluster connected to the electrodes through me- 
chanically soft organic linkers. A periodic self-oscillation 
of the cluster in conjunction with the cluster charg- 
ing/decharging is predicted for sufficiently large bias volt- 
ages, leading to an average current proportional to the 
self-oscillation frequency. 3 Two factors neglected in the 
original proposal may complicate the analysis and pre- 
vent the observation of the shuttle effect in molecular- 
assembled single-electron devices^ (1) At small oscilla- 
tion amplitude, the cluster displacement modulates the 
metal-molecule bond length/strength rather than the 
molecule core. The exponential dependence of the tun- 
neling resistance on the displacement of the metal island 
may not hold; 7 (2) At large oscillation amplitude, the in- 
ternal structure of the linker molecule may be distorted 
due to the forces induced by the change in metal-molecule 
bond, the net charges on the cluster and the applied elec- 
trical field. The shuttle mechanism of electron transfer 
has also been used to interpret recent experiment on a 
Cqo single-electron transistor,-' 9 but alternative explana- 
tions^ also exist. Recently it has been suggested that 
measurement of shot noise spectrumAi and full count- 
ing statistics^ may help in elucidating the mechanism 
involved. 

A simpler way of incorporating mechanical degree of 
freedom into single-electron devices is to couple one of 



the tunnel junctions to a nanomechanical oscillator ji 3 . us- 
ing e.g., the microcantilever tip of a conducting Atomic 
Force Microscope (AFM) or hybrid Scanning Tunneling 
Microscope (STM)/AFMi 14 i 15 i 16 The model system is il- 
lustrated schematically in Fig. ^a) , where the nanopar- 
ticle is connected to the substrate through an organic 
molecule (assumed rigid), while tunneling across the top 
contact can be modulated mechanically through the vi- 
brating tip. The introduction of a vibrating tunnel con- 
tact leads to tunneling and displacement currents flow- 
ing simultaneously in accordance with the probe vibra- 
tionfl£>l& which can be measured separately using a two- 
phase lock-in amplifier^ The discrete-electron tunneling 
also induces an electrostatic force on the tipiiii^iiSiiS 
Since both the tunneling/displacement currents and the 
force on the microcantilever tip can be measured as a 
function of the bias voltage, useful information regard- 
ing the interplay among tip vibration, discrete electron 
motion and metal-molecule interaction can be extracted, 
which may potentially allow applications in displacement 
detection and chemical/bio-sensing. The purpose of this 
paper is thus to present a theoretical analysis of such 
nanomechanical double-barrier tunneling junctions using 
realistic atomic-scale models. 

Master- equation approach to single- electron tunneling 
through molecular-assembled metallic nanoparticles. — 
We consider periodic vibration of the tip, with the tip- 
nanoparticle distance being x(t) = do + d\cos2iT ft. The 
coupled molcculc-nanoparticle- vibrating tip system is de- 
scribed by the following Hamiltonian: 



H = H S + H T + H! + V es 
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where H a = Y,k e k a a\ a a ka {a = S,T) and Hi = 

eiicjjCu describe the noninteracting electrons in the 
substrate (S), tip (T) electrodes and on the central island 
(I) respectively. We model both the electrodes and the 
nanoparticle as infinite electron reservoirs with electro- 
chemical potential fi a ,a — S,T,I. The indices k and I 
enumerate the electron states of the electrodes and the 
nanoparticle. We describe the linker molecule using an ef- 
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FIG. 1: (Color online) (a) shows schematic illustration of 
the nanomechanical double-barrier tunnel junctions formed 
using a vibrating platinum tip and a gold metallic nanoparti- 
cle connected to the gold substrate through an alkane dithiol 
molecular bridge with 12 alkane units (AK12). The plat- 
inum tip is modeled as a 4-atom pyramid sitting on top of the 
semi-infinite substrate, (b) shows the equivalent circuit model 
of the system. Note that the substrate junction capacitance 
(Cs) and resistance (Rs) are fixed, while the tip junction ca- 
pacitance (Ct) and resistance (Rt) are time-dependent due 
to the tip vibration. The tip-substrate capacitance (Co) is 
also time-dependent. 



fective single-particle Hamiltonian H mo i = Yli e ibl m bim- 
The electrostatic part of the energy is V es (qs,qT) — 

2 

5t^+qsVs+qTV T , C s = C S +C T , where C S {t) and q S ( T ) 
are the capacitance and charge of the substrate (tip) 
junction respectively^ The voltage drop across the sub- 
strate (tip) junction is determined from the tip-substrate 

bias voltage V by Vs(t) = ~gr^ V- The net charge on 
the central island is q — qs — qr = — e(n + n x ), where 
en x is the background charge and en is the quantized is- 



land charge. Unlike the conventional single-electron de- 
vices where the background charge is typically induced 
by charged impurities embedded in the insulating layer, 
here an intrinsic background charge n x can be induced by 
the charge transfer between the molecule and the central 
island during the formation of the substrate junction, 21 
and can be obtained from microscopic molecular junction 
calculations^ 

The transfer of single-electrons is mediated by tunnel- 
ing through the molecular (substrate) junction Ht s = 
J2ik( t i^iksbl m a k s+U m -,iibl m ciie^' t ' s +H.C.), and the tip 

junction H TT (t) = 12ik( t ll-M x ( t ))"kT c lT e ~ i * T + H.C.), 
where we have shown explicitly the time-dependence 
of the tip-nanoparticle coupling through the tip posi- 
tion x(t). The phase operator <frs(T), canonically con- 
jugate to the charge qscr) 011 the substrate (tip) junc- 
tion, keeps track of the quantized electron tunneling since 
&S(T),qs(T} = ie and e^ s ^q s(T} e-^ s m = qs(T) - 



FIG. 2: (a) shows the molecular junction and tip junction 
resistance as a function of the junction length. The tip junc- 
tion resistance increases exponentially with tip-nanoparticle 
distance as the tip remains weakly coupled to the nanoparticle 
during the tip vibration cycle, (b) shows the static current- 
voltage (I-V) and excess island charge-voltage (n-V) charac- 
teristics of the double-barrier tunnel junctions for a 10(nm)- 
diameter nanoparticle at tip-nanoparticle distance of 3.0A 
(dotted line), G.0A (solid line) and 9.0 A (dashed line) for tem- 
peratures of T = 5, 150K respectively. The inset shows the 
magnified view of I-V characteristics at 6.0A and 9. OA Note 
that at x = 6. OA, we already have Rt » Rs, so the n-V 
characteristics remains approximately the same as x increases 
further. 
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FIG. 3: (a) shows the time-averaged tunneling I-V char- 
acteristics (solid line) and the instantaneous tunneling I-V 
characteristics at x — 3.0A (dotted line) at three tip vibrat- 
ing frequencies of / = 1, 10, 100MHz. The tunnel currents 
are virtually independent of the tip vibrating frequency, (b) 
shows the time-averaged n-V characteristics (solid line), the 
instantaneous n-V characteristics at x = 3.0j4 (dotted line) 
and x — 9.0A (dashed line) respectively. The instantaneous 
values are the values obtained at the moments when the tip 
moves to x — 3.0A or 9.0A 
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FIG. 4: Upper figures show the average number of shut- 
tled electron per cycle of tip vibration, middle figures show 
the average displacement current, and lower figures show the 
average electrostatic force on the metallic tip as a function 
of tip-substrate bias voltage at tip vibrating frequencies of 
/ = 1, 10, WO(MHz). In the middle figures, we have shown 
both the capacitive part (solid line) and "Coulomb Blockade" 
part (dotted line) of the displacement current. The capac- 
itive part increases linearly with voltage and tip vibrating 
frequency. But the "Coulomb Blockade" part oscillates with 
increasing voltage. In the lower figures, we have shown both 
the capacitive part Foap (solid line) and discrete-electron part 
Fcb (dotted line) of the electrostatic force on the tip. Note 
that the step- wise increase in the shuttled electrons in the up- 
per figures correlates with the step-like increase in discrete- 
electron force Fcb in the lower figures. 



£^24 rp-^g gi ec t ros tatic energy change after tunneling of 
an electron through the substrate(tip) junction is thus 

A n = e^s(T)T/ ese - l s( T) _ Ves = eVfr ~ — i— is^LL 2 

q 2 20.23 



'S(T) 
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We describe single-electron transport using the master 
equationi 4 i 18 i 25 



dP m (t) 
dt 



P m +i (*)r~ (A m , t) + p m _ x (t)r+ ( A m _! , t) 

P m (t) (T-(A m _ 1 ,t)+ T+ ( A m , t) ) (2) 



where P m is the probability of finding m excess elec- 
trons on the nanoparticle and r + (~' = 



-+(-) 



The rate for transition m — ► mil due 
to electron tunneling through the substrate (tip) junc- 
tion, r^y, is obtained from the electron trans- 
mission coefficient of the corresponding junctions us- 
ing T+ {T) (w,t) = I J dET siT) (E,V;t)f s(T) (E)(l - 

fi(E-w)) and T- {T) (w,t) = \ J dET s(T) (E, V ; t) (1 - 
fs(T){E))fi(E - w), 21 where w = A m or A m _i. The 
transmission coefficient T$(t) can be determined from the 
standard Non-Equilibrium Green's Function theory2& as 
developed for molecular tunnel junctions. 21-22 Here the 



time dependence in the tunneling rate T + ( ~> is intro- 
duced through the time-dependent tunneling rate across 
the tip junction 1^ ' due to the tip vibration x{t). Since 
the tunneling time for electron tunneling across the tip 
junction is much smaller than the tip vibration period, 
the time-dependence in Tr indicates that we calculate the 
transmission coefficient using the tunneling Hamiltonian 
HT T (t) at the instantaneous tip position x(t). 

The tunneling current across the tip junction and the 
average excess electron number on the nanoparticle at 
time t are given by I^ un (t) = e£ m P m (t)[T+(A m ,t) - 
r~(A m _!,i)] and n(t) = J2 m mP m(t) respectively. The 
displacement current flowing through the vibrating tip is 
due to time-variation of system charges and is fixed by 
electrostatic consideration a 18 ' 23 : 



I° lS (t) = C V+^= I CaP (t) + I CB {t) 



Cs + Ct 



)V + ei 



Ct 



dVCs + Cr 



(n(t) + r43) 



where Co is the capacitance between the tip and the sub- 
strate. The first term leap is the displacement current 
induced by the direct capacitive coupling between the 
tip and the substrate and the series capacitance of the 
double-barrier junction. The second term Iqb gives the 
"Coulomb blockade" component of displacement current 
due to the time variation of the excess electrons on the 
nanoparticle. The electrostatic force on the oscillating 
tip can also be divided into the capacitive force Fc a p 
and the discrete-electron force Fcb as follows: 



Fi 



tip 



lrfCp 2 l d(q T V T ) 
2 dx 



1 d 



2 dx 

-c T (gm 2 



Fcap + 

e d 
~ 2dx { 



Fcb 
CsCt 



(n + n x )lW) 



Since our main interests here are in electron transport 
and charging- induced forces, we neglect non- voltage de- 
pendent forces like the long-range van der Waals force. 
Since we consider transport within the Coulomb blockade 
regime, the tip remains weakly coupled to the nanopar- 
ticle during the entire cycle of vibration, and we neglect 
also the short-range contact force that may arise from 
the tip-nanoparticle bonding. Both type of forces would 
have been important if we want to study the detailed 
oscillation dynamics of the cantilever^ 

Device model. — We apply the theory to single-electron 
tunneling through a 10(nm)-diameter gold nanoparti- 
cle connected to the gold (111) substrate through an 
alkane dithiol molecule with 12 alkane units (AK12) and 
probed using a vibrating platinum tip. Electron trans- 
port through gold-dithiol molecule-gold junctions have 
been recently studied in detail using a first-principles self- 
consistent matrix Green's function theory^ from which 
we obtain the intrinsic background charge en x induced 
by the charge transfer from the molecule to the nanopar- 
ticle and the transmission coefficient for electron tun- 
neling through the metal- molecule-nanoparticle junction. 
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To summarize, the gold electrodes are modeled as semi- 
infinite (111) single-crystals. Six nearest-neighbor gold 
atoms on each metal surface (twelve gold atoms overall) 
are included into the "extended molecule" where the self- 
consistent calculation is performed. The rest of the elec- 
trodes (with the six atoms on each surface removed) are 
considered as infinite electron reservoirs, whose effects 
are included as self-energy operators. The calculation is 
performed using the Becke-Perdew-Wang parameteriza- 
tion of density-functional theory— and appropriate pseu- 
dopotentials with corresponding optimized Gaussian ba- 
sis setsj2£ We find that the resistance of the AK12 molec- 
ular junction is Rs ~ 300(MQ), and the background 
charge is n x « — 0.45^ 

Since the tip geometry in STM/AFM is often not 
well characterized, we model the Pt tip as a 4-atom Pt 
pyramid sitting on top of the (111) Pt substrate. The 
tunneling matrix elements are obtained using the semi- 
empirical Extended Huckel Theory (EWT)— and consid- 
ering coupling between the Pt pyramid and six neighbor 
gold atoms on the surface of the gold nanoparticle, with 
the apex Pt atom sitting in front of the center of the 
six gold atoms arranged following the (111) gold surface 
geometry (Fig. ^a)). The calculated tip junction re- 
sistance Rt as a function of the tip apex-nanoparticle 
distance is shown in Fig. HJa), which increases expo- 
nentially with the tip-nanoparticle distance once it falls 
within the weak-coupling regime. For capacitance mod- 
eling, such atomic-scale model is not needed. The tip 
is instead replaced by a 1/xra-radius sphere, which takes 
into account approximately the tip curvature effect and 
is common in AFM research. 16 - 30 Here both the substrate 
and tip junction capacitance are obtained from classical 
electrostatics of conducting sphere sitting in front of the 
conducting substrate^ which gives the substrate junc- 
tion capacitance of 0.9aF, the tip junction capacitance 
of 0.97aF and tip-substrate capatitance of 0.34/i* 1 at tip- 
nanoparticle distance x = 9.0A respectively. 

Results and their interpretation. — We assume the tip 
vibration being described by oscillation amplitude of 
di = 3.0 A and average tip-nanoparticle distance of do = 
6. OA The device functions as a nanomechanical single- 
electron turnstile, which can be understood by examin- 
ing the static (tip not moving) excess island electron- 
voltage (n-V) characteristics at different tip-nanoparticle 
distances x (Fig. 0(b)). We find that the trend of the 
charging state on the nanoparticle is reversed as the tip- 
nanoparticle x increases from 3.0A to 9.0A. This is be- 
cause at x = 3.0A, the tip junction resistance is much 
smaller than the substrate junction resistance while at 
x = 9.0A, the tip junction resistance becomes much 
larger than the substrate junction resistance (Fig. 13(a)). 
As the tip starts vibrating at x = do + d\ = 9.0A, the 
nanoparticle will be charged by discrete-electron tunnel- 
ing across the substrate junction if the bias voltage is 
large enough to overcome the charging energy. As the tip 
moves close to the nanoparticle at x — dg — di = 3.0A, the 
tip junction becomes sufficiently conductive (Rt << Rs) 



allowing the stored excess charge on the nanoparticle to 
be transferred onto the tip. A discrete number of elec- 
trons is thus being shuttled across the system per tip 
vibrating cycle. Note that the n-V characteristics at 
x = do = 6.0A is similar to that at x = 9.0A because 
Rt is already much larger than Rs at x = 6.0A, which is 
also the case during most part of the tip vibration cycle. 

We solve the rate equation (Eq.[2J numerically, 31 from 
which we obtain the time-dependent tunneling current, 
the displacement current and the electrostatic force on 
the Pt tip. We calculate the probability P/\(m) that 
m electrons have been transferred during one cycle of 
tip vibration from Pa (to) = ^2kPi{k)P(k — m\k,to + 
1//) (f is the vibrating frequency)^ where Pi(k) is the 
possibility of finding k electrons on the nanoparticle at 
time to- P(k—m\k, to+1/ /) is the conditional probability 
of finding k — to electrons on the nanoparticle after one 
cycle of tip oscillation if there are k electrons at initial 
time to, which is obtained by solving the rate equation 
Eq. [3 with the initial condition P(k,to) — l,P(i,to) = 
0, i kX The average number of electrons being shuttled 
across the system is n shu ttie = J2 m mP A (m). 

The average (time-average over one cycle of the tip 
vibration) tunneling current-voltage (I-V) flowing across 
the tip junction and excess island electron-voltage (n-V) 
characteristics of the vibrating double-barrier junctions 
at temperature of 5K are shown in Fig. [3] for three tip 
vibrating frequencies f — 1, 10, lQO(MHz). For compar- 
ison, we have also shown the instantaneous tunneling I-V 
characteristics when the tip is closest to the nanoparticle 
(x = 3.0A) in Fig. OJa), the instantaneous n-V char- 
acteristics at two different tip-nanoparticle distances of 
x = 3.0, and 9.0A in Fig. |3fb). The instantaneous val- 
ues at x = 3.0A, 9.0A are the values obtained at the 
moment when the tip moves to x = 3.0A, 9.0A respec- 
tively. Since the instantaneous tunneling current of the 
vibrating double-barrier junction decreases rapidly as tip 
moves away from the nanoparticle, and the tip velocity 
approaches zero as tip moves closer to the nanoparti- 
cle, the average tunneling I-V characteristics is similar in 
both magnitude and voltage dependence to the instan- 
taneous tunneling I-V characteristics at x — 3.0 A at all 
three tip frequencies. Note that the tunneling currents 
shown here are virtually independent of the tip vibrating 
frequency, since the time it takes for electron to tunnel 
across the tip junction is the smallest of all time scales in- 
volved here and can be considered as instantaneous com- 
pared to the tip vibration. 

The n-V characteristics show stronger tip-frequency 
dependence than the tunneling I-V characteristics. The 
instantaneous n-V characteristics at x = 3.0A changes 
little with tip frequency. They are also similar to the 
static n-V characteristics in Fig. [2Jb). This is because 
at x=3 A, tip junction resistance is small. The corre- 
sponding RC charging time is much smaller than the tip 
vibrating period, the instantaneous charging state on the 
nanoparticle is thus not affected by the tip vibration. As 
the tip moves away from the nanoparticle, Rt increases 
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rapidly leading to larger RC charging time. Therefore 
charge imbalance on the nanoparticle can be affected at 
large tip-nanoparticle distance as the tip vibrating fre- 
quency increases, which is seen clearly from the instanta- 
neous n-V characteristics at the largest tip-nanoparticle 
distance x — 9. OA The time-averaged n-V characteris- 
tics approximates well the average between the instanta- 
neous n-V characteristics at x = 3.0A and 9.0A 

Fig. 0] shows the number of electrons n s huttie being 
shuttled across the system per cycle of tip vibration, the 
displacement current and the electrostatic force on the 
tip averaged over one cycle of tip vibration at tip vibrat- 
ing frequencies of / = 1, 10, 100MHz respectively. Note 
that n a huttie decreases with increasing tip frequency due 
to reduced efficiency of charge transfer from the central 
island to the tip. The capacitive part of the displacement 
current leap increases linearly with the tip-substrate bias 
voltage (solid line), but the "Coulomb Blockade" part 
of the displacement current Iqb oscillates with the bias 
voltage, since dn/dt changes sign during one cycle of the 
tip vibration. The peak and dip positions in Iqb corre- 
late closely with the step- wise change in the shuttled elec- 
trons with respect to the bias voltage. The magnitude of 



leap increases linearly with increasing tip vibration fre- 
quency, but the increase in magnitude of Iq b is less reg- 
ular. We have also shown the time-averaged capacitive 
force Fcap and the discrete-electron force Fqb on the tip 
separately in Fig. 0] Note that Fc a p is quadratic voltage- 
dependent and is attractive. In contrast, the discrete- 
electron force Fqb is repulsive. The magnitude of Fcb 
shows step-wise increase with bias voltage as more elec- 
trons are allowed to be shuttled, similar to the voltage- 
dependence of n shuttle*^ The overall magnitude of Fcap 
and Fcb is comparable. Unlike the displacement cur- 
rent, the magnitude of the discrete-electron force depends 
only weakly on the tip frequency through the frequency 
dependence of n shu ttie- 

To conclude, we have presented a microscopic the- 
ory of mechanical modulation of single-electron tunnel- 
ing through a molecular-assembled metal nanoparticle in- 
duced by a vibrating nanoprobe. The effects discussed 
here may be of interest for applications in nanomechani- 
cal chemical/bio-sensors. 

This work was supported by the DARPA Moletron- 
ics program, the DoD-DURINT program and the NSF 
Nanotechnology Initiative. 
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